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^— «j ' Abstract. 

A reduction of total mean turbulent pressure due to the presence of magnetic fields was pre- 
viously shown to be a measurable effect in direct numerical simulations. However, in the studied 
^ parameter regime the formation of large-scale structures, as anticipated from earlier mean-field 
simulations, was not found. An analysis of the relevant mean-field parameter dependency and 
the parameter domain of interest is conducted in order to clarify this apparent discrepancy. 
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' Strong magnetic fields at the solar surface are generally thought to originate from the 

C/^ . coherent rise of magnetic flux tubes from the tachocline through the solar convection 

^ ' zone. While the idea is elegant, the question remains whether it can be considered as 

Q_i, more than a toy model, as the physics of the creation and rise of these flux tubes is 

o 
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not sufliciently understood ( Parkeill2009f ). Thus far any 'successful' numerical simulation 



of this process had to rely on strong assumptions, be it in th e initia l conditions or in 
simplified equations such as the thin fiux tube approximation ( Spruitl[l981i) . As such it 



, makes sense to explore alternative mechanisms of magnetic structure formation. Several 
models have been proposed where large -scale magnetic field c oncentra tions are created 



, through instabilities at the solar surface (iKitchatinov fc Mazur 2000; Br andenburg et al 

> ■ E)10aL hereafter BKR). 

i Turbulence is generally associated with enhanced transport effects. However, it can 

' also generate structures on much larger scales than its drivin g scales; see as an exam 



' pie the inverse cascade in 2D hydr odynamic turbule nce (jKraichnan 1967 ) or 3D MHD 



turbulence with magnetic helicity (jPrisch et al. 1975 ). We study here the interaction of 



, the turbulence with a background magnetic field. From the approximate conservation of 
' total turbulent energy j^tot j BKR find a reversed feedback from the magnetic fluctuations 



on the turbulent pressure (jRogachevskii fc Kleeorinll2007l hereafter RK): 



^tot- 



It can be seen that the effective mean magnetic pressure force is reduced and can be 
^ . reversed in a certain parameter range. It was suggested by RK that this positive feedback 
" " ' could lead to an instability, resulting in the concentration of magnetic flux. 

BKR confirmed the validity of approximate turbulent energy conservation using di- 
rect numerical simulations (DNS) of homogeneous isothermal turbulence and they also 
demonstrated the basic phenomenon of magnetic flux concentration through the inter- 
action between turbulence and the mean Lorentz force in mean-field MHD simulations 
which led to a linear instability for sufficiently strong stratification. These results were 
recently corroborated using DNS of inhomogeneous (stratified) isothermal turbulence of 
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Figure 1. Left: Normalized effective mean magnetic pressure as a function of deptli for 
Bo = O.lBeqo, Bp = . 2i?eao, and Bo — O.SBcqo using Re = 120 and Vtm = 1- Adapted 
from [Brandenburg et alj (|2010bf ). Right: Visualization of By — Bp on the pe r ipher y of the com- 
putational domain Bo = O.lBoqO, and Ptm ~ 2. Adapted from lKemel et aLl ([20101 ). 



a cubic computational domain. An example is shown in Figure [T] where we show that the 
normalized effective mean magnetic pressure does indeed become negative in a large part 
of the domain for different imposed field strengths. However, the formation of large-scale 
structures is not observed, as can be seen from the right-hand panel of Figure [TJ 

It is important to understand whether there is really a conflict between DNS and 
mean-field results, or if these two studies simply apply to different parameter regimes. 
In order to have a chance to resolve this question, we conduct a systematic parameter 
survey of the instability in the mean-field model. More specifically, we determine the 
functional dependence of the growth rate on the input variables in order to find the 
relevant parameter space for the instability to develop. 

In an isothermal stratified box we solve the equations in two dimensions: 

V In p + g + J^M + ^K,tot, 
qp)W] + B-V[{l-qs)B] 

^K,tot = (J^t + 1^) (V^tJ + VV • 17 + 2SV Inp) 

is the total (turbulent and microscopic) viscous force with S being the viscous stress 
tensor. In addition, we solve the continuity and uncurled induction equations, 

^ = -[/. Vlnp- V-C7, ^ = C7 X B - (ryt + r?) J. 

We adopt a Cartesian coordinate system, (x, y, z). The mean field is given by B = 
(0, Bq, 0) + V X a, and the vertical gravitational acceleration is g — (0, 0, —g). The other 
input parameters of the simulations are the sound speed Cs, the density at the top of the 
box ptop, the pressure coefficient qp (RK), the magnetic Prandtl number PrM and the 
molecular diffusivity rj. 

We start by isolating the dependence of the growth rate on the turbulent diffusivity 
and find approximately linear behaviour. From a data survey we observe, depending on 
the applied field, different diffusivity values above which the instability does not develop. 



where 

p^M = -5V[(l 
is the mean-field Lorentz force and 
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Figure 2. Left: dependence of growth rate on turbulent diffusivity and viscosity. Right: 
modified growth rate as a function of $. 

which would imply 

A = - (77t, + z/t)fc' 

where 1/fc is a length scale introduced for dimensional reasons. We find that $ is propor- 
tional to the Alfven speed at the top of the box, where the instability initiates. Eventually 
we arrive at the fit formula 

A + (l + Pr'y'^) Tjtk^ = WAfc, (1 + qp/q;) exp (4/i?p) , 



as shown in Figure H Here, Pr^j,"^'*^ 



vt/vt, is some typical vertical length scale, and 



Hp ~ c^/g is the density scale height. The instability resulting from this feedback effect 
was verified but not observed to generate large-scale structures in DNS with the current 
scale separation and parameter range studied so far. 

An extension of this work will be the inclusion of more physics (e.g. radiative cooling) 
in the DNS and the comparison with mean-field models in a 3D setup. As was already 
pointed out in BKR, the mean-field model predicts additional short-wavelength pertur- 
butions along the direction of the mean magnetic field, that were not included in the 
present study. 
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